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This paper discusses particle production in Schwarzchild-like spacetimes and in an uniform electric 
field. Both problems are approached using the method of complex path analysis which is used to 
describe tunnelling processes in semi-classical quantum mechanics. 

Particle production in Schwarzchild-like spacetimes with a horizon is obtained here by a new and 
simple semi-classical method based on the method of complex paths. Hawking radiation is obtained 
in the (t, r) co-ordinate system of the standard Schwarzchild metric without requiring the Kruskal 
extension. The co-ordinate singularity present at the horizon manifests itself as a singularity in the 
expression for the semi-classical propagator for a scalar field. We give a prescription whereby this 
singularity is regularised with Hawking's result being recovered. The equation satisfied by a scalar 
field is also reduced to solving a one-dimensional effective Schrodinger equation with a potential 
( — l/x 2 ) near the horizon. Constructing the action for a fictitious non-relativistic particle moving in 
ON ' this potential and applying the above mentioned prescription, one again recovers Hawking radiation. 

GN . In the case of the electric field, standard quantum field theoretic methods can be used to obtain 

particle production in a purely time-dependent gauge. In a purely space-dependent gauge, however, 
the tunnelling interpretation has to be resorted to in order to recover the previous result. We 
attempt, in this paper, to provide a tunnelling description using the formal method of complex 
f"*^ ' paths for both the time and space dependent gauges. The usefulness of such a common description 

becomes evident when 'mixed' gauges, which are functions of both space and time variables, are 
analysed. We report, in this paper, certain mixed gauges which have the interesting property that 
mode functions in these gauges are found to be a combination of elementary functions unlike the 
standard modes which are transcendental parabolic cylinder functions. 
QQ ' Finally, we present an attempt to interpret particle production by the electric field as a tunnelling 

£SJ ' process between the two sectors of the Rindler spacetime. 
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I. INTRODUCTION AND SUMMARY 
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, In this paper, we present a critical review of particle production in Schwarzchild-like spacetimes and in an uniform 
electric field in Minkowski spacetime. We approach both problems by the method of complex paths discussed by 
Landau in |^| where it is used to describe tunnelling processes in non-relativistic semi-classical quantum mechanics. 
This powerful technique will be used as the basis to provide a new method of recovering Hawking radiation in the 
usual Schwarzchild co-ordinates without requiring the Kruskal extension. It will also be used to describe particle 
production in an electric field in different non-trivial gauges and to link particle production in an electric field to 
tunnelling processes occuring in the Rindler spacetime. 

The Schwarzchild-like spacetimes we consider here are the usual black hole, the Rindler and the de-Sitter spacetimes. 
In the standard black-hole spacetime, particle production was obtained by Hartle and Hawking Q] using semi-classical 
analysis. In this method, the semi-classical propagator for a scalar field propagating in the Schwarzchild spacetime is 
analytically continued in the time variable t to complex values. This analytic continuation gives the result that the 
probability of emission of particles from the past horizon is not the same as the probability of absorption into the 
future horizon. The ratio between these probabilities is of the form 

P[emission] = P[absorption]e _ ' 3f ' (1-1) 

where E is the energy of the particles and f3 = l/87rM is the standard Hawking temperature. The above relation 
is interpreted to be equivalent to a thermal distribution of particles in analogy with that observed in any system 
interacting with black body radiation. In the latter case, the probability of emission of radiation by the system is 
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related to the probability of absorption by the system by a similar relation as given above. In Hawking's derivation, 
the Kruskal extension is of vital importance in obtaining the thermal spectrum. 

In this paper, we propose an alternate derivation of Hawking radiation without using the Kruskal extension. Our 
motivation for using the standard Schwarzchild co-ordinates rather than the Kruskal system are the following: (1) 
The Schwarzchild spacetime is a static spacetime. It contains a global Killing vector £ a and the symmetry generated 
by this vector is respected in the Schwarzchild co-ordinates (t, r). That is, the surfaces t — constant where t is the 
Schwarzchild 'time' variable (this Killing vector is space-like in the region interior to the horizon but it is still a 
symmetry of the system) has the same structure for all t; g a b is independent of t. (2) The surface area of spheres of 
constant 'radial' co-ordinate r happen to be Airr 2 , which is that of a sphere in flat spacetime, and hence these surfaces 
can be used to measure r. In contrast, the Kruskal system is an explicitly time- dependent spacetime since metric 
components now depend on the Kruskal time co-ordinate. It does not possess a global Killing vector associated with 
surfaces u = constant where u is the Kruskal 'time' co-ordinate. Indeed, the presence and structure of such a global 
Killing vector is not immediately apparent from the form of the metric in these co-ordinates. Surfaces of constant 
'radial' co-ordinate do not possess the same property that corresponding surfaces in the Schwarzchild spacetime do. 
Further, the sectors which contain the past horizon and the time-reversed copy of the usual Schwarzchild spacetime 
are unphysical. In any realistic collapse scenario, these sectors cannot exist. The non-static nature of the Kruskal 
manifold imply that the 'in' and 'out' vacua cannot be the same and explicit particle production can take place. 

The difficulty, of course, is that the standard Schwarzchild co-ordinates possess a co-ordinate singularity at the 
horizon. As we shall see, this bad behaviour of the co-ordinates appears as a singularity in the expression for the 
semi-classical propagator near the horizon and we have to provide a specific prescription to bypass it. This prescription 
gives the same result as that obtained by Hawking and can be used in all spacetimes with a Schwarzchild-like metric. 
Note that the method described above is fundamentally different from the usual method of calculating the Bogolubov 
coefficients for an eternal black hole given, for example, in Ref. 0. In Ref. Q], appropriate linear combinations of the 
Schwarzchild mode functions that are analytic on the full Kruskal manifold (except at the past and future singularites 
at r = 0) are constructed. The scalar field is then expanded in terms of these modes with the vacuum state being 
the Kruskal vacuum. Such an expansion provides the appropriate connnection between the annihilation and creation 
operators for the scalar field in the Schwarzchild and Kruskal co-ordinate systems. Then, for a given time like observer 
in either the left or right Schwarzchild sector, whose vacuum state is the Schwarzchild vacuum, the number operator is 
easily calculated and is found to correspond to a thermal spectrum. In our method, however, the action functional is 
constructed using the Hamilton- Jacobi method in the appropriate co-ordinates. The singularity in the action caused 
by the singular behaviour of the co-ordinate system (in the unextended Schwarzchild co-ordinates) at the horizon is 
handled by the prescription to obtain particle production. We also reduce the problem of the massive scalar field 
propagating in such a spacetime to an effective Schrodinger equation which has a singular effective potential near 
the horizon. In this case, the singular behaviour of the co-ordinates manifests itself as a singularity in the potential. 
The non-relativistic action for a particle moving under the influence of such a potential is constructed and the same 
prescription is used to bypass this singularity. We again recover Hawking radiation. 

We next consider the problem of a scalar field propagating in flat spacetime in an uniform electric field background. 
The total particle production due to the presence of the electric field upto the one-loop approximation is correctly 
calculated by the gauge invariant method proposed by Schwinger jlj. The same problem can be reduced, in a time 
dependent gauge, to an equivalent Schrodinger equation with an inverted harmonic oscillator potential. Such an 
equation can be solved exactly using the standard flat spacetime quantum field theoretic methods. Since the problem 
is explicitly time dependent, the vacuum state at t — ► — oo and at t — > oo are not the same. The Bogoliubov coefficients 
between the "in" and "out" vacua are easily calculated and the total particle production turns out to be the same as 
that calculated by the Schwinger method. However, if a space dependent gauge is used to describe the same field, the 
vacuum state of the field remains the same for all time and hence no particle production can take place. To recover 
the standard result, it is usual to resort to the tunnelling interpretation. This interpretation is useful since it provides 
a dynamical picture of particle production (see, for example 0) and is the only way in which the standard gauge 
invariant result can be recovered in a time independent gauge. 

In this paper, we attempt a tunnelling description based on the method of complex paths for both the time-dependent 
and time-independent gauges. This method is used to calculate the transmission and reflection coefficients (or the 
tunnelling coefficients) for the equivalent quantum mechanical problem. Then, an interpretation of these coefficients, 
in order to explicitly obtain the standard gauge invariant result is provided. However, the usefulness of the tunnelling 
description is seen when simple 'mixed' gauges, which are functions of both space and time, are considered. In some 
gauges, the scalar wave equation can be reduced to solving effective Schrodinger equations in suitable new variables 
which are combinations of the usual spacetime co-ordinates. By applying the method of complex paths to the these 
Schrodinger equations, the tunnelling coefficients are seen to match those obtained either for the purely time or 
purely space dependent gauges. The gauge invariant result is now recovered using the appropriate interpretation 
needed to identify the Bogolubov coefficients. In certain other gauges, the mathematics becomes simpler and the field 
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equation can be reduced to a first order differential equation rather than an effective Schrodinger equation cannot be 
made. For such cases, the exact mode functions themselves are used to set up a tunnelling scenario and the resulting 
tunnelling coefficients are interpreted according to the tunnelling interpretation to recover the gauge invariant result. 
The mixed gauge functions we report in this paper have the interesting and useful property that the mode functions 
are combinations of elementary functions. This is in contrast to the standard modes which are transcendental due to 
the presence of the parabolic cylinder functions. These new modes are found to be singular on the lightcone. This 
property is very similar to the modes of the Schwarzchild-like spacctimes since these too are singular on the horizon. 

In classical theory, the action of an uniform electric field on a charge imparts a constant uniform acceleration to 
it. The spacetime metric in the rest frame of the charge is the Rindler frame. Quantum field theory on the other 
hand, predicts particle production arising due to the presence of an electric field in the spacetime. It is therefore of 
interest to ask if particle production is linked in some way to the the presence of a Rindler frame since both are, in 
a sense, natural for this problem. We attempt to link particle production by an uniform electric field with processes 
occuring in the Rindler frame by proposing an interpretation of the standard result in terms of tunnelling between 
the two Rindler sectors. We do this in a heuristic manner and show that this tunnelling process between th e R indler 
sectors gives rise to the leading exponential factors in the expression for the effective lagrangian (see Eqn. ( [f.l|) ) ||. 

The layout of the paper is as follows. Section (||) contains the semi-classical derivation of Hawking radiation 
without taking recourse to the Kruskal extension. It also contains the reduction to an effective Schrodinger equation 
and the subsequent recovery of Hawking radiation in the limit of large black hole mass. Section (III) discusses particle 
production in an uniform electric field in a set of new and non-trivial gauges. It relies heavily on the method of 
complex paths outlined in Appendix (A). In section (IV), we attempt to link particle production in the electric field 
to processes occuring in the Rindler frame. 



II. PARTICLE PRODUCTION IN SPACETIMES WITH HORIZON 



Hawking's result that a black hole radiates is essentially a semi-classical result with the thermal radiation arising 
because of the presence of a horizon in the spacetime structure. We will review briefly the conventional derivation 
of the thermal radiation using path integrals. Consider a patch of spacetime, which in a suitable co-ordinate system, 
has one of the following forms (we assume c = 1): 

ds 2 = B{r)dt 2 - B-^rjdr 2 - r 2 {d6 2 + sin 2 (9) d<t> 2 ) (2.1) 



ds 2 = B{x)dt 2 - B- 1 {x)dx 2 - dy 2 - dz 2 (2.2) 

where B(r) and B{x) are functions of r and x respectively. The horizon in the above spacetimes is indicated by the 
surface r = ro (x ~ xq) at which B(r) (B(x)) vanishes. We furthe r as sume that B'(r) ~ dB/dr (B'(x) = dB/dx) is 



finite and non-zero at the horizon. Co-ordinate systems of the form (2.1) can be introduced in parts of the Schwarzchild 
and de-Sitter spacetimes while that of the form 



flat spacetime. Given the co-ordinate system of ( 



(2.2) with the choice B{x) = 1 + 2gx represents a Rindler frame in 
(2.1) say, in some region R, we first verify that there is no physical 
singularity at the horizon, which in the case of the Schwarzchild black hole, is at the co-ordinate value ro = 2M where 
M is the mass of the black hole. Having done that, we extend the geodesies into the past and future and arrive at 
two further regions of the manifold not originally covered by the co-ordinate system in ( |2.1| ) (the Kruskal extension). 
It is now possible to show that the probability for a particle with energy E to be lost from the region R in relation 
to the probability for a particle with energy E to be gained by the region R is given by the relation 

P. — P p~ l3E (1 3 s ! 

J loss — ± gain c V / 

where f3 — 8irM. This is equivalent to assuming that the region R is bathed in radiation at temperature f3~ . In the 
derivation given in the paper by Hartle and Hawking Q , thermal radiation is derived using the semiclassical kernel by 
an analytic continuation in the time co-ordinate t to complex values and it was shown that the probability of emission 
(loss) from the past horizon was related to absorption (gain) into the future horizon by the relation (j^) . 

Since all the physics is contained in the (t, r) plane (or the (t, x) plane), we will discuss Hawking radiation in 1 + 1 
dimensions first and show subsequently that the results generalise naturally to 3 + 1 dimensions without modification. 
We first derive the semi-classical result in the (t,r) (or {t,x)) plane by applying a certain prescription to bypass the 
singularity encountered at the horizon. After this, we reduce the problem of the Klein-Gordon field propagating in 
the Schwarzchild spacetime to an effective Schrodinger problem in (1+1) dimensions and rederive the semi-classical 
result by using the same prescription. 
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A. Hawking Radiation in 1+1 dimensions 



Consider a certain patch of spacetime in (1+1) dimensions which in a suitable co-ordinate system has the line 
element, (with c = 1) 



B(r)dt 2 - B^(r)dr 2 



(2.4) 



where B(r) is an arbitrary function of r. We assume that the function B(r) vanishes at some r — r$ with B'(r) = 
dB/dr being finite and nonzero at r$. The point r — ro indicates the presence of a horizon. It can be easily verified 
that no physical singularity exists at the horizon since the curvature invariants do not have a singularity on the 
horizon. Therefore, near the horizon, we expand B(r) as 



B{r) = B'{r )(r - r ) + 0[(r - r ) 2 ] = R{r Q )(r - r ). 



(2.5) 



where it is assumed that R(ro) + 1 0. We now use the equation satisfied by the minimally coupled scalar field $ 
with mass mo propagating in the spacetime represented by the metric ( [2.4] ) to obtain the Hamilton- Jacobi equation 
satisfied by the action functional Sq. The semi-classical propagator can be constructed using So which will be used to 
analyse the singularity at the horizon. (We emphasise again that this method is different from that used to compute 
the Bogolubov coefficients using appropriate superpositions of mode functions as outlined in 0.) 

The equation satisfied by the scalar field is, 



where the □ operator is to be evaluated using metric (2.4). Expanding the LHS of equation ( |2.6| ), one obtains 



1 <9 2 $ d 



B(r) dt 2 dr 



Or 



(2.6) 



(2.7) 



The semiclassical wave functions satisfying the above equation are obtained by making the standard ansatz for 
which is, 



$(r, t) = cxp 



S(r,t) 



(2.8) 



where S is a function which will be expanded in powers of h. Substituting into the wave equation (2/7), we obtain, 



B(r) V dt 



OS 



os 



dr 



1 d 2 S _ £PS_ _ dB{r) OS 
B~{r)W~ dr~~dr 



(2.9) 



Expanding S in a power series of {fi/i), 



S(r,t) = S (r,t) + - SiM) + - S a (r,t) 



and substituting into Eqn fl2.9|) and neglecting terms of order (h/i) and greater, we find to the lowest order, 



(2.10) 



B(r) 



dSo 
dt 



B(r) 



dSo 
dr 







(2-11) 



Eqn (2.11) is just the Hamilton- Jacobi equation satisfied by a particle of mass mo moving in the spacetime determined 
by the metric (2.4). The solution to the above equation is 



S (r,t) 



-Et± 



dr 



E 2 - m 2 B{r) 



(2.12) 
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where E is a constant and is identified with the energy. Notice that in the case of mo = 0, Eqn (2.11) can be exactly 
solved with the solution 



S (r,t;m = 0) = F 1 (t-r*)+F 2 (t + r*) (2.13) 
where the "tortoise" co-ordinate r* is defined by 

r ^!wr (2 - 14) 

and Fi and F 2 are arbitrary functions. If Fx is chosen to be F\ = —Et + Er* and F 2 chosen to be F 2 = —Et — Er*, 



then it is clear that the solution given in Eqn (2.13) is the same as that in Eqn (2.12) with mo set to zero. Therefore, 
in the case mo = 0, the semiclassical ansatz is exact. In the following analysis we will specialise to the case mo = 
for simplicity. The case mo ^ will be considered later. The essential results do not change in any way. 

The semiclassical kernel K (r 2 , t 2 ; n, t\) for the particle to propagate from (£i,ri) to (t 2 ,r 2 ) in the saddle point 
approximation can be written down immediately as. 

K{r 2 ,t 2 ;n,ti) =iVexp (jrS {r 2 ,t 2 -,ri,ti)\ (2.15) 

where So is the action functional satisfying the classical Hamilton- Jacobi equation in the massless limit and N is a 
suitable normalization constant. So(r 2 , t 2 ; r\, t\) is given by the relation 



r T2 dr 

S (r 2 ,t 2 ; n, ti) = Sb(2, 1) = ~E(t 2 - h) ± E — - 

Jri B ( r ) 



(2.16) 



The sign ambiguity (of the square root) is related to the "outgoing" (dSo/dr > 0) or "ingoing" (dSo/dr < 0) nature 
of the particle. As long as points 1 and 2, between which the transition amplitude is calculated, are on the same side 
of the horizon (i.e. both are in the region r > ro or in the region r < ro), the integral in the action is well defined 
and real. But if the points are located on opposite sides of the horizon then the integral does not exist due to the 
divergence of B~ l (r) at r = rg. 

Therefore, in order to obtain the probability amplitude for crossing the horizon we have to give an extra prescription 
for evaluating the integral || . Since the point B = is null we may carry out the calculation in Euclidean space or 

cquivalcntly — use an appropriate it prescription to specify the complex contour over which the integral has to be 
performed around r — r . The prescription we use is that we should take the contour for defining the integral to be an 
infinitesimal semicircle above the pole at r = ro for outgoing particles on the left of the horizon and ingoing particles 
on the right. Similarly for ingoing particles on the left and outgoing particles on the right of the horizon (which 
corresponds to a time reversed situation of the previous cases) the contour should be an an infinitismal semicircle 
below the pole at r — ro. Equivalently, this amounts to pushing the singularity at r — ro to r — ro =F * e where the 
upper sign should be chosen for outgoing particles on the left and ingoing particles on the right while the lower sign 
should be chosen for ingoing particles on the left and outgoing particles on the right. For the Schwarzchild case, this 
amounts to adding an imaginary part to the mass since ro = 2M. 

Consider therefore, an outgoing particle (dSo/dr > 0) at r = r\ < ro. The modulus square of the amplitude for 
this particle to cross the horizon gives the probability of emission of the particle. The contribution to So in the ranges 
( r ij r o — e) and (ro + e, r 2 ) is real. Therefore, choosing the contour to lie in the upper complex plane, 

r +e 



f r ° +t dr 

So [emission] = —E lim / + (real part) 



iirE , _. 
= -^r^ + real part 2.17 

where the minus sign in front of the integral corresponds t o th e initial condition that dSo/dr > at r = n < ro 



For the sake of definiteness we have assumed R(ro) in Eqn (2.5) to be positive, so that B(r) < when r < ro. (For 
the case when R < 0, the answer has to be modified by a sign change.) The same result is obtained when an ingoing 
particle (dSo/dr < 0) is considered at r = r\ < r . The contour for this case must be chosen to lie in the lower 
complex plane. The amplitude for this particle to cross the horizon is the same as that of the outgoing particle due 
to the time reversal invariance symmetry obeyed by the system. 

Consider next, an ingoing particle (dSo/dr < 0) at r = r 2 > ro. The modulus square of the amplitude for this 
particle to cross the horizon gives the probability of absorption of the particle into the horizon. Choosing the contour 
to lie in the upper complex plane, we get, 
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f r °~ e dr 

Sq [absorption] = —E lim / — — + (real part) 
^0/ B(r) 



ittE , , 

= — + (real part) (2.18) 

-R(ro) 

The same result is obtained when an outgoing particle (dSo/dr > 0) is considered at r = T2 > rg. The contour for 
this case should be in the lower complex plane and the amplitude for this particle to cross the horizon is the same as 
that of the ingoing particle due to time reversal invariance. 

Taking the modulus square to obtain the probability P, we get, 

f 2ttP\ 

P[emissionJ oc exp — — — - I (2-19) 



and 



so that 



hR 



P[absorption] cx exp ( ] (2.20) 



/ 4nE\ r . 
P[emissionJ = exp — - — I P [absorption]. (2-21) 

\ hR J 

Now time reversal invariance implies that the probability for the emission process is the same as that for the absorption 
process proceeding backwards in time and vice versa. Therefore we must interpret the above result as saying that the 
probability of emission of particles is not the same as the probability of absorption of particles. In other words, if the 
horizon emits particles at some time with a certain emission probability, the probability of absorption of particles at 
the same time is different from the emission probability. This result shows that it is more likely for a particular region 
to gain particles than lose them. Further, the exponential dependence on the energy allows one to give a 'thermal' 
interpretation to this result. In a system with a temperature the absorption and emission probabilities are related 
by 

P[emission] = exp(— /3P)P[absorption] (2.22) 



Comparing Eqn (2.22) and Eqn (2.21), we identify the temperature of the horizon in terms of R(ro). Eqn (2.21) is 
based on the assumption that R > 0. If R < there will be a change of sign in the equation. Incorporating both the 
cases, the general formula for the horizon temperature is 



/T 1 = M (2.23) 

47T 



For the Schwarzchild black hole, 



/ 9M\ 1 

B{r) = (l - — J « — (r - 2M) + 0[(r - 2Mf] (2.24) 

giving R = (2M) _1 , and the temperature /3 _1 = fr/8irM. For the de-Sitter spacetime, 

B(r) = (1 - PV) w 2H(H- 1 - r) = -2H(r - P _1 ) (2.25) 

giving P^ 1 = TiH /2tt. Similarly for the Rindler spacetime 

B(r) = (1 + 2gr) = 2g(r + (2g)~ 1 ) (2.26) 

giving (3~ l = gh/2ir. The formula for the temperature can be used for more complicated metrics as well and gives 
the same results as obtained by more detailed methods. 

The prescription given for handling the singularity is analogous to the analytic continuation in time proposed by 
Hawking Q] to derive Black hole radiance. If one started out on the left of the horizon and went around the singularity 
r = ro by a 2n rotation instead of a rotation by tt, it can be easily shown that it has the effect of taking the Kruskal 
co-ordinates (v,u) to (—v,—u). A full rotation by 2ir around the singularity can be split up into two parts to give 
the amplitude for emission and subsequent absorption of a particle with energy E. Since the amplitudes for the two 
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processes are not the same in the presence of a horizon, one obtains the usual Hawking radiation given in Eqn (2.21) 
with the value of R(r a ) being (2M)" 1 . This process is similar to that given in 14] which relates the amplitudes 
involving the past and future horizons. In Hawking's paper, analytically continuing the time variable t to t — AMiir 
takes the Kruskal co-ordinates (v, u) to (— v, —u) and since the path integral kernel is analytic in a strip of AMiir 
below the real t axis, Hawking radiation is obtained by deforming the contour of integration appropriately. Note 
however — in our approach — we did not require the Kruskal extension and worked entirely in the (t, r) co-ordinates. 

When mo_^_0, the validity of the semi-classical ansatz must be verified. To do this, consider the perturbative 
expansion ( [2.1C| ). Retaining the terms of order Ti/i and neglecting higher order terms, one finds, upon substituting 
for Sq given by the relation (2.12) and solving for Si, 



Si = -Eit ± EEi 



B(r) y/E 2 - m 2 B(r) 4 v y " 



(2.27) 



where E\ is a constant. From the above equation, it is seen that S\ has a singularity of the same order as So at r = Tq. 
When calculating the amplitude to cross the horizon, the contribution from the singular term just appears as a phase 
factor multiplying the semiclassical kernel and is inconsequential. The non-singular finite terms do contribute to the 
kernel but they contribute the same amount to S [emission] and S [absorption] and they do not affect the relation 
between the probabilities P [emission] and P [absorption]. Subsequent calculation of the terms 5*2, S3, and so on, show 
that all these terms have a singularity at the horizon of the same order as that of So- Their contribution to the 
probability amplitude is just a set of terms multiplied by powers of h which can be neglected. From this we can 
conclude that the semiclassical ansatz, in the perturbative limit, is a valid one. 

The generalization to (3+l)-dimensions is straigh tforw ard. We will work with (2.1) which is in spherical polar 



co-ordinates. The results obt ained are extendable to (2.2) in a straightforward manner. The Klein-Gordon equation, 
written using the metric (2.1), is 



1 2 $ 1 / 2 „. ,0$ 
3«Z r B ( r )«T 



B(r) dt 2 r 2 dr 



1 



dr ) r 2 sin(6>) 89 

1 2 $ 



d ( . d$ 
sm(60 — 



r 2 sm z (0)d^ % z 

Since the problem is a spherically symmetric one, one can put $ = \&(r, t)Y l m (6, <fi) to obtain, 



1(1 + 1) m 2 



1 d 2 ^ 19/ „„. . 0* 

r B(r) 

B{r) dt 2 r 2 dr\ K ' dr 

Making the ansatz W = exp((i/h)S(r,t)) and substituting into the above equation, we obtain, 



(2.28) 



(2.29) 



1 



B(r) 

h 
+- 

1 



dS_ 
~dt 

1 d 2 S 



B{r) 



dS_ 

dr 



i{i + i)h 2 



B ^9 2 S 1 d{r 2 B) dS 



= 



B(r) dt 2 dr 2 r 2 dr dr 

Expanding S in a power series as in Eqn ( [2.10| ), we obtain, to the zeroth order in h/i, 

dr 



So 



-Et± 



B(r) 



^E 2 -B{r)(m 2 +L 2 /r 2 ) 



(2.30) 



(2.31) 



where L 2 = 1(1 + l)h 2 is the angular momentum. It is easy to see from the above equation that near the horizon, 
the pres ence of the L 2 term can be neglected since it is multiplied by B(r). Therefore, the semiclassical result of 
section (II A) follows even in the case of (3+l)-dimensions. The semi-classical ansatz is valid in this case as can be 
seen by calculating explictly the higher order terms in the expansion for S. All these terms have a singularity at the 
horizon of the same order as that of So and they contribute to the semiclassical propagator either as phase factors or 
as terms multiplied by powers of h~ which are entirely negligible. Expanding the Klien-Gordon equation for $ using 
the metric (2.2) gives analogous results and will not be explictly given. 
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B. Reduction to an effective Schrodinger Problem in (1 + 1) dimensions 



Consider the relativistic equation for the wave function $ in Eqn (2.7). We include the mass mo here but we shall 
see later that it does not appear in the final result. Setting 



*(r) 



-iEt/H 



Q(r) 



we get the equation, 



d 2 Q(r) 
dr 2 



B"{r) (B'(r)f 



E 2 



2B{r) AB 2 (r) h 2 B 2 (r) h 2 B(r) 



Q(r) = 



(2.32) 



(2.33) 



where B' = (dB/dr) and B" = (d 2 B/dr 2 ). Near the horizon, we use the expansion of B{r) given in Eqn ([2.5]). 
Neglecting terms of order l/(r — ro) as compared to terms of order l/(r — ro) 2 , we get, in the limit of fi — > 0, 



d 2 Q(r) 
dr 2 



(r - r ) 2 



Q(r) = 



where 



E 2 



fJ 



h 2 R 2 



(2.34) 



(Notice that mo does not appear to the leading order in the above equation. Very close to the horizon, the term 
containing the mass does not contribute significantly. Eqn (2.34) is therefore applicable to both massless and massive 
scalar particles.) Making the transformation x = (r — ro), we finally obtain the effective Schrodinger equation for the 
system with h = 2m = 1 with a potential (—g/x 2 ). 



d 2 Q(x) 



4qw = 



(2.35) 



dx 2 x" 

This potential is symmetric but singular at the origin x = 0. To make the analogy with the schrodinger equation, we 
will replace the righthand side of Eqn. (2.35) by EQ(x) and finally take the limit of E — > 0. This "energy" E should 
not be confused with the energy E of the field in the original relativistic system. We, therefore, consider 



d 2 Q{x) 
dx 2 



X 



EQ(x). 



(2.36) 



The energy spectrum is continuous for all values of E which, for E < 0, is peculiar to this potential since, for energies 
less than the potential energy, the spectrum is usually discrete. 

The semiclassical analysis follows closely the method adopted in section (II A). The action functional A for a 
classical particle moving in a potential —g/x 2 satisfies the Hamilton- Jacobi equation 



dA 
~dt 



dA 

dx 







The solution can be immediately written down as, 

A = -Et ± 



— \/£'< 2 

X 



(2.37) 



(2.38) 



Eqn (2.38) has an integral which is divergent if the action is computed for points lying on the opposite sides of the 
horizon x = 0. Since this has a similar form to Eqn ( 2.12| ), the prescription used in evaluating Remission] and 
5 [absorption] can be similarly used to evaluate A [emission] and .4 [absorption]. The results are 



^[emission] = in^/g + (real part) 
^[absorption] = — i^^fg + (real part). 



(2.39) 



Constructing the semiclassical propagator as before and taking the modulus square to obtain the probabilities for 
outgoing and ingoing particles, we get 



P [outgoing] = exp 



4nE 

~hR 



P [ingoing] 



(2.40) 
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The temperature /3 _1 for the system is the same as that in Eqn ( 2.23Q and one recovers the usual result. 

To verify that the semi-classical analysis is valid, one must compute the correction terms and check that these have 
a singularity of the same order as possessed by A. To do this, consider the effective Schrodinger equation ( 2.36 ) with 
factors of Ti put in. 



%2 d 2 Q(x) 
dx 2 



4<2(*) = EQ{x) 



Putting Q(x) — exp(iA(x)/?i), and substituting into Eqn (|2.41 ), 



-ifi 



d 2 A(x) 
dx 2 



( dA{x) 
V dx 



= E 



Expanding A in powers of h/i, we get, 



A = A+ -A 1 



-) A 2 

i 



(2.41) 



(2.42) 



(2.43) 



Substituting into Eqn ( 2.42| ) and proceeding as usual, we find that A is given by Eqn ( 2.38| ). The next term A\ is 
given by 



M =g 



dx 



1 



x Ex 2 + g 



(2.44) 



The relation for Ai also has a singularity at the origin of the same order as A. Explicit calculation of the subsequent 
terms in the expansion of A reveals that all these terms have a singularity of the same order as that of A and therefore 
their net contribution to the kernel is either as phase factors or as the exponential of finite terms multiplied by powers 
of h. Therefore, we conclude as before that the semiclassical ansatz is valid. 

We show now that the effective Schrodinger equation in (3+l)-dimensions is the same as in Eqn ( 2.35 ). We 
conside r he re the reduction of the Klein-Gordon equation in spherical polar co-ordinates obtained in Eqn (2.28) using 
metric (2T). Setting 



and substituting into Eqn ( [2.28 ), we obtain, 



d 2 * 1 d(r 2 B) d^ 



E 2 m 1 L 



dr 2 ' r 2 dr dr ' \h 2 B{r) Ti 2 ti 2 r 2 



* = 



(2.45) 



(2.46) 



Making the substitution 



we get the result, 



d 2 Q 
dr 2 



1 

B 2 



(B 



i\2 



E 2 



y/r 2 B{r) 



Q(r) 



1 

B 



B' 



L 2 
h 2 r 2 



= 



(2.47) 



(2. 



where B' and B" are th e fir st and second derivatives of B{r) respectively. Near the horizon r = ro, using the expansion 
for B(r) given in Eqn (2.5), it is easy to see that the 1/B 2 term in the above equation dominat es ove r the 1/B term 
which can therefore be neglected. The resulting Schrodinger equation is the same as in Eqn ( 2.35 ) in the limit of 
h — > 0. It can easily be proved that the effective Schrodinger equation for the cartesian metric (2.2) gives exactly the 
same result. 







III. PARTICLE PRODUCTION IN AN UNIFORM ELECTRIC FIELD 



We now move onto the discussion of particle production in an uniform electric field. We study a system consisting 
of a minimally coupled scalar field $ propagating in flat spacetime in an uniform electric field background. The two 
standard gauges, namely the purely time dependent and the purely space dependent gauges, are first considered and it 
is shown how the method of complex paths outlined in Appendix (Ala) can be used to construct a viable tunnelling 
interpretation in each case. The standard quantum field theoretic result will not be rederived here since it is well 
known (though its results are used to construct the tunnelling interpretation). The tunnelling description is then 
applied to a few non-trivial but simple 'mixed' gauges which are functions of both space and time. 

The method of complex paths is a useful tool used to calculate the transmission and reflection coefficients in 
semi-classical quantum mechanics in one space dimension || and we briefly summarise it here, leaving the details 
to Appendix (A). (Readers unfamiliar with this approach should read Appendix (A) at this juncture since those 
results will be used extensively). First two disjoint regions are identified where the semi-classical wave functions 
can be written down. (The Schrodinger potential that will be encountered most often will be the inverted harmonic 
oscillator potential (— x 2 ) which has the semi-classical, disjoint, regions located at x — > ±00.) Then a tunnelling 
scenario is set up by imposing appropriate boundary conditions on the semi-classical wave functions in these regions. 
One region is assumed to contain the transmitted wave while the other contains the incident and reflected waves. To 
obtain the tunnelling coefficients, the solution is analytically continued to the complex plane and the behaviour of the 
transmitted wave is studied along a complex contour (with the space variable now considered complex) joining the 
two regions. The contour is chosen such that the semi-classical condition is satisfied all along the contour. Rotation 
along the contour transforms the transmitted wave either to the incident wave or to the reflected wave thus relating 
the transmission amplitude to either the incident or the reflected wave amplitude. Using the normalisation condition 
satisfied by the tunnelling coefficients, both coefficients can be determined. 

The complex contour should be chosen so that singularities (real and complex) in the potential, where the semi- 
classical ansatz is invalid, are avoided. If such singularites exist, they contribute to the determination of the tunnelling 
coefficients. Notice that this method works for the exact mode functions too (as it must). An appropriate tunnelling 
scenario has to be set up and the complex path method can be applied. Such situations will be encountered below 
in certain gauges where the explicit modes are known but the problem cannot be reduced to an effective Schrodinger 
equation. 



A. Time dependent gauge 

The four vector potential giving rise to a constant electric field in the x direction is assumed to be of the form 

A" = (0, -E t, 0, 0) (3.1) 

The electric field is E = E$k. The minimally coupled scalar field $ propagating in flat spacetime, satisfies the 
Klein-Gordon equation, 



((0 M + iqA^){d" + iqA»)+m 2 )§ 







(3.2) 



where m is the mass and q is the charge of the field. The mode functions of $ can be expressed in the form 
$(t, x) = /k(i)e lk ' x where fu(t) satisfies the equation, 



dt 2 

Introducing the variables, 



we obtain the equation, 



/ k + [m 2 + k\ + (k x + qE t) 2 ] / k = ; kj. = {k y , k z ). 



<qE t 



f qEa 



k\ + m 2 
qE Q 



(3.3) 



(3.4) 



- r 2 / k = A/ k 



(3.5) 



The above equation is essentially a Schrodinger equation in an inverted oscillator potential with a positive "energy" 
A. Since the e nergy is positive, the problem is essentially an over the barrier reflection problem. Using the results of 
section (Alb), we can calculate the reflection and transmission coefficients exactly as 
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where we have put h — 2m = g\ = 1 and set E\ = A in Eqns. (A24, A25). To identify the Bogoliubov coefficients a\ 



and /3\, we recast the normalisation condition R + T = 1 in the form, 

± - £ = 1 (3-7) 
and then identify \f3\\ 2 with R/T and |q;a| 2 with l/T. Therefore, the Bogoliubov coefficients are given by, 

|^| 2 = e --=expf-^±^ 
V qE 

K| 2 = e-- + 1 = exp(-^l±^l) + 1. (3.8) 

The transmission and reflection coefficients are time reversal invariant and are dependent only on the energy (magni- 
tude and sign). They are also independent of the direction in which the boundary conditions are applied. To obtain 
a dynamical picture of particle production, we have to interpret these quantities suitably. In the present case, the 
following interpretation seems adequate. A purely positive frequency wave with amplitude square T in the infinite 
past, t — ► — oo, evolves into a combination of positive and negative frequency waves in the infinite future t — ► oo with 
the negative frequency waves having an amplitude square R and the positive frequency waves having an amplitude 
unity. The quantity R/T determines the overlap between the negative frequency modes in the distant future and 
the positive frequency modes in the distant past (the notation here differs from the treatment given in || , Q ) . This 
is identified with the modulus square of the Bogoliubov coefficient f3\ which is the particle production per mode A. 
Using the normalisation condition satisfied by the Bogoliubov coefficients, |a,\| 2 — \(3\\ 2 = 1, |o;a| 2 can be calculated 
to be l/T. Once the Bogoliubov coefficients have been identified, the effective lagrangian can be easily calculated. 
This derivation will not be repeated here. We refer the reader to Ref. j|] and Ref. || for the explicit calculation. 
Note that the particular interpretation given in this case is due to its similarity with the more rigourous calculation 
by quantum field theoretic methods. In the next section, in which we discuss the space dependent gauge, we will be 
forced to adopt a different interpretation in order to identify particle production. 



B. Space dependent gauge 



The four vector potential giving rise to a constant electric field in the x direction is now assumed to be of the 
form 



A>* = (-E x, 0, 0, 0) 



(3.9) 



The electri c fie ld is E = E gx a s before. The field <& satisfies Eqn. ([T^) as before. Substituting for the potential A^ 
from Eqn. (3J)) into Eqn. ( |3.2| ), we obtain, 

(3.10) 
(3.11) 



(<9 2 - V 2 - 2iqE xd t - q 2 E 2 x 2 + m 2 ) $ = 

We write $ in the form 

$ = e - lult e tk y y+lk * z (f>(x) 
and obtain the differential equation satisfied by cj> as 

^ + ((w + qE x) 2 -kl-m 2 )cj) = 
where we have used the notation k\ = k 2 + k 2 . Making the following change of variables 



'qE x 



/ qEo 



k\ + m 2 

qE 



(3.12) 



(3.13) 



into the differential equation for (f>, it reduces to the form, 



11 



"ay 



p 2 (f> = -\<j> 



(3.14) 



In this form, we see that the above differential equation has the form of an effective Schrodinger equation with an 
inverted harmonic oscillator potential and an negative energy — A. If we apply the results of section (Alt), we obtain 
the result for tunnelling through the barrier. Following the treatment in Ref. pj and using the results of section (Alb) 
we can calculate the reflection and transmission coefficients exactly as 



R = 



T = 



(3.15) 



where we have put h = 
R + T = 1 in the form 



2m = gi = 1 and set E\ = —A in Eqns. (A24, A25). We cast the renormalisation condition 



1 

R 



T 
R 



(3.16) 



and then identify the rate of particle production per mode with T/R. The interpretation of particle production using 
the tunnelling picture now proceeds as follows. A right moving travelling wave of amplitude square 1/ R is incident on 
the potential. A fraction T/R is transmitted through it and a wave of unit amplitude is scattered back. The tunnelling 
probability, which is T/R, is interpreted as the rate at which particles are being produced by the background electric 
field. This matches exactly with the expression for \f3\\ 2 given in Eqn. (3.8). With this interpretation, we recover the 
usual gauge independent result. 



C. Mixed gauges 

We shall now study the problem in a new set of gauges which prove to be useful and instructive. (As far as the 
authors know these gauges have not been studied in the literature before.) When parameters which specify these 
gauges are varied, the problem can be mapped either onto a "tunnelling through the barrier" or "over the barrier 
reflection" system. Then, using the tunnelling description developed previously for each of these systems, the gauge 
invariant result is obtained. 

For certain ranges of these parameters, the scalar wave equation can be reduced to solving a second order equation 
which can be convert ed to an effective Schrod inger equation. This equation is studied using the complex path 



method as in section ( [II A ) and section ( [IIB ). The solutions to these effective Schrodinger equations are usually 
transcendental in form. 

For other parameter ranges, however, the scalar field equation reduces to a first order differential equation whose 
solution is a combination of elementary functions. To recover the gauge invariant result, the solution and its complex 
conjugate are used to set up a tunnelling scenario from which the tunnelling coefficients are obtained. Then, using 
the tunnelling interpretation, the gauge invariant result of particle production is recovered. 



1. Gauge Type 1 

The first gauge type we consider is a simple generalisation of the space-dependent gauge (in Eqn. (|]^)) of the form 

A M = (-E x + Dit, 0, 0, 0) (3.17) 

where Eq is as usual the magnitude of the electric field and D\ is a gauge parameter which can be positive or negative. 
It is easily verified that the above gauge type gives the same electric field as the pure time and space gauges did. The 
differential equation for the scalar field $ is, 

(d 2 - V 2 + 2iq(-E x + D^dt - q 2 (-E x + Drf) 2 +iqD 1 +m 2 )$ = (3.18) 

We make a judicious choice of variables of the form 

u = -E Q x + Dit ; v = E Q x + D\t. (3.19) 



Notice that the coefficients in Eqn. ( 3.1 8| ) are dependent only on the variable u and not on v, y or z and so the above 
equation is seperable in the variables (u,v,y, z). Expressing the derivatives (dt,d x ) in terms of (d u ,d v ) and writing 
4> in the form 
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(3.20) 



one obtains 



(Df El) 



du 2 



2i(qD lU -j(Di + E$)) 
+ (m 2 + k\ + iqD x - 



d(j> 
du 



El 



2 2 

q u 



2qDau) (j>{u) = 



(3.21) 



It is easy to see from the above equation that two distinct cases can be identified here, namely, \D\\ ^ \Eq\ and 
D\ = ±Eq. In the first case, the differential equation is a second order one and the effective Schrodinger equation 
can be obtained by eliminating the first derivative. In the second case, however, the resulting equation is a first order 
differential equation whose solution is an elementary function. 
Consider first the case \Di\ ^ |-Eo|- Writing (f>(u) in the form, 



Q(u) exp- 



qD x 



1{D\ 



El) 



2{D{-E 2 ) 



El) 



and defining a new variable p - 



, — (2jDi/q), one obtains the effective Schrodinger equation as, 
d 2 Q q 2 E 2 p 2 ^ _ (m 2 + k 2 ± ) 



dp 2 



(D 2 -E 2 ) 2 



Q(p) 



(Df 



El 



Q(p)- 



(3.22) 



(3.23) 



The effective potential is clearly seen to be that of an inverted oscillator. Notice that, depending on whether D\ < El 
or D\ > El, the problem reduces to a "tunnelling thro ugh the barrier" probl em o r an "over the barrier reflection" 
problem respectively. Using the results of section (Alb) and especially Eqn. (A21), the value of the quantity 2e\ is 
found to be (with U = 2m = 1) 



2ei = 



-E x = 
9i 



qE 



-Sgn(Dl-El) = \Sgn(L%-E*) 



(3.24) 



where Sgn(x) is the sign function which is positive for x > and negati ve fo r x < and A is defined in Eqn. ( |3.5| ). 
Thus, the tunnelling coefficients are given by either Eqn. (3.6) or Eqn. ( 3.1 5| ) depending on Sgn(D 2 — El) but are 
neverthless independent of D\ as expected. The tunnelling interpretation required to recover the standard result 
proceeds accordingly. Note that only the Q(u) part of the full solution to $ contributes to the transmission and 
reflection coefficients. The solutions to Q(u) are the usual parabolic cylinder functions 

Now, consider the more interesting case D\ = 
reduces to 



-Eq with u = Eo(t — x). The differential equation Eqn. (3.21) for 



2i(qE u-2 1 El)^- + (m 2 + k\ + iqE - q 2 u 2 + 2qE ol u) <j>{u) 
du 







The solution is easily obtained to be 

4>{u) - 




i\/2 - 1/2 



-iqu /4E 



(3.25) 



(3.26) 



with A being defined in Eqn. (3.13). This solution resembles the asymptotic forms for some of the parabolic cylinder 
functions except that it is exact and is clearly a combination of elementary functions. (Notice that the solution is 
singular on the surface t — x = 2"//q which is remniscent of the behaviour of black hole modes near the horizon. The 
implications of this will be discussed in a future publication.) To recover the gauge invariant result in this case, notice 
that the comple x conj ugate <jf (u) is also a solution. With this pair of independent modes, one can apply the theory 



given in section ( A 1 a| ) to set up a tunnelling problem with the appropriate boundary conditions at u = ±oo. This is 



most conveniently done by defining a new dimensionless variable 



(3.27) 



The mode function d> now becomes 



13 



») = ( S -2 7v ^7?) 



iA/2 - 1/2 



-is 2 /4 



(3.28) 



with an analogous expression for 4>*(s). We now assume that the wave function is a right moving travelling wave 4>r 
as s — ► +00 while as s — > —00, it is a superposition of an incident wave of unit amplitude and a reflected wave given 
by (f>L. Therefore, we have, 



t>R = Cga-^-VSe" 3 /* 



= (_ S )'V2-l/2 e -i S 2 /4 + C2( _ s) -a/2-l/2 e i s ; 



/4 



+00) 
— 00) 



(3.29) 



where C3 and C2 are the transmission and reflection amplitudes respectively. We have made the approximation 
s — 2 r fi l VEo/q w s which holds in the limit s —> ±00. Applying the method of complex paths (see Appendix (A)) to 
rotate (f>R about the upper complex contour (in accordance with the semi-classical ansatz) , one obtains 



C*2 



-iC'3 exp 



ttA 
~2~ 



Using the normalisation condition | C 2 1 2 + | C3 1 2 = 1, we get, 



R 



|C 2 | 2 



1 



T 



1 



(3.30) 



(3.31) 



which are the usual tunnelling coefficients given in Eqn.( 3.15 ). It is therefore seen that the system c orresp onds to an 
tunnelling through the barrier problem. The interpretation, of course, follows that given in section ([II B). 
For the case Z?i = — Eq, the modes are given by 



<P(u) 




iA/2 - 1/2 



Jqu 2 /4E 



(3.32) 



An analysis similar to the case D\ = +Eq shows that the tunnelling coefficients are the same as in Eqn. ( |3.3l| ) Thi s 
system also corresponds to tunnelling through the barrier with the corresponding interpretation given in section (III B) 
required in order to recover paticle production. Therefore, it is seen that, for the case \D±\ 7^ \Eq\, the magnitude of 
D\ decides the appropriate tunnelling interpretation even though it does not appear in the final expressions for the 
tunnelling coefficients. In contrast, for the cases D\ = ±Eq, the system reduces to a tunnelling through the barrier 
with the sign of D\ not playing any role in deciding the appropriate interpretation. 

The above considerations also work for an analogous generalisation of the time-dependent gauge of the form 



= (0,-E o t + D lX , 0, 0, 0) 
in an obvious manner and we will not repeat the discussion. 



(3.33) 



2. Gauge Type 2 

The second gauge type we consider is of the form, 

= {Dit - D 2 x, D 2 x — D\t, 0, 0) (3.34) 

where D\ and D 2 are arbitrary constants such that E = Di + D 2 - The magnitude and direction of the uniform 
electric field is seen to be the same. It is in this gauge that, for D\ = D 2 , the mode functions have the simplest form 
possible. Writing out the differential equation for the scalar field $ and defining new variables 

u = D 1 t-D 2 x ; v = D 1 t + D 2 x (3.35) 

and setting $ = e lk yV+ lk ^ z e -n" <fi(u) as before, one obtains, 

{Dl-Dl)^± + K[ q (D 1+ D 2 )u- 7 (Dl+Dl)] d ± 

+ [to 2 + k 2 ± + iq{D 1 + D 2 ) - 7 2 (D 2 - D 2 ) + 2 1 q{D 1 - D 2 )u] 4>{u) = (3.36) 



14 



Here too, two distinct cases can be distinguished, namely, D\ ^ D2 and Di = D2 (we will discuss the case D2 
which corresponds to a 'pure gauge' with zero electric field a little later). 
Consider first the case D\ ^ £>2- Writing 4> in the form, 



exp - 



n 2 ) 



2{Di-D 2 ) (£>? 

and introducing a new variable p = u — (2"/DiD2/qE ), one obtains an effective Schrodinger equation as, 

„2 7^2 J2 



d 2 Q 
dp 2 



(D 2 -D 2 ) 2 



(3.37) 



(3.38) 



where we have used Eq = D\ + £>2- This equation has the same form as that in Eqn. (3.23). Depending on whether 
D\ < D 2 or D\ > £)§, the problem reduces to a "tunnellin g thro ugh the barrier" prob lem or an "over the barrier 
reflection" problem respectively. Using the results of section (Alb) and especially Eqn. (A21), one finds the value of 
the quantity 2ei to be (with h = 2m = 1) 



2e, = 



1 k\ + n 
Ei = 

gi qE Q 



-Sgn{D\-Dl) = \Sgn(Df- D\) 



(3.39) 



Thus, the tunnelling coefficients are given by either Eqn. (3.6) or Eqn. (3.15) depending on Sgn(D 2 — D%) and are 
dependent only on Eq as expected. The tunnelling interpretation required to recover the standard result proceeds 
accordingly. 

However, when D\ = D2, the differential equation for 4>(u) reduces to 



Ai{qD 1 u- 1 D\) < ^ + (m 2 + k 2 ± + 2iqDi)4>(u) 







with the solution 



<j>{u) 




i\/2 - 1/2 



(3.40) 



(3.41) 



This is the simplest mod e function possible for any gauge of the electric field. This solution is also singular just like 
the modes in Eqn. (3.26) but on the surface t — x = 7/(7. The form of the function <P{u) clearly indicates the facto r 
responsible for the non-zero tunnelling coefficients. One can, in a manner similar to that outlined for the modes (3.26), 
use the independent mode functions <fi and <jf to set up a tunnelling problem in the limit u — > ±00. Rotating the 
transmitted wave in the upper complex plane (this is the only co ntou r possible in accordance with the semi-classical 
ansatz) the tunnelling coefficients are found to be those in Eqn. (3J3). The system, in this case, correspo nds to over 

the barrier reflection and the interpretation required to obtain particle production follows that in section ( HI A y 

Finally, we look at the case D2 — —D\. The electric field, in this pure gauge, is identically zero. Solving Eqn. (|3.36|) 
for this case, the solution for 6 is found to be 



4>{u) 



_ -iXu/2-fD-L -iqu 2 /2D 1 



(3.42) 



It is clear from the above form of <f> that the transmission coefficient is unity and the reflection coefficient is zero. No 
particle production takes place. Explicit calculations, of course, verify this. The tunnelling interpretation gives a null 
result in this case. 

The mixed gauge types in Eqn. ( 3.17 ) and Eqn. ( ]3 .34j) represent the same electric field. Hence, these gauges must 
be related to the pure space and time gauges in Eqns ( |3.l| , 3.9 ) by gauge transformations of the form A a — A a + d a f 
where / is a suitable gauge function. Choosing the base gauge A a as the space-dependent gauge in Eqn. ( |3.9| ), one 
obtains the following gauge functions for the two mixed gauge types. 



A a — ► A a = {-E x + Dit, 0, 0, 0) / = - D Y t 2 

A a -> A a = {Dit - D 2 x, D 2 x - Dit, 0, 0) / = D x tx - 



Dit 2 



Dox 2 ) 



(3.43) 



The gauge transformed scalar fields i n th e A a gauge ought to be relat ed to those in the A a gauge by a phase factor 



But the solutions give in Eqn. ( 3.23j ) (with Di = Eq) and Eqn. ( 3.41 ) (with D\ = D% = Eq/2) are clearly not 
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gauge transformations of the mode functions of gauge (3.9) which contain parabolic cylinder functions. These modes 
are intrinsic to the gauges themselves and cannot be obtained by a simple gauge transformation of the modes of the 
space-dependent gauge. This is reminiscent of the situation in the Rindler and Minkowski frames. The Rindler plane 
wave modes of the form e~~ lWT /(£) (where r and £ are the time and space co-ordinates in the Rindler frame) are not 
obtained from the Minkowski plane wave modes of the form e ~ luJt + kx (where t and x are the Minkowski time and 
space co-ordinates) by a co-or dinate transform ation . 

The solutions given in Eqn. (|]|| and Eqn. can be obtained in a more general and elegant fashion as follows. 

Consider the Minkowski metric expressed in the (u, v, y, z) co-ordinate system where u = t — x and v = t + x. 



ds 2 



dudv 



dy 2 



dz A 



Assume a vector potential of the form 



A» = (g(u), h(u), 0, 0) 



(3.44) 



(3.45) 



where g(u) and h(u) are arbitrary functions with A 1 - 1 = (A u , A v , A V ,A Z ) being expressed in the (u, v, y, z) co-ordinate 
system and not in the (t, x,y, z) system. The electric field is given by E = (dg(u) / <9it)x. Notice that the v component 
of the vector potential, A v — h(u), does not contribute to the electric field and is a 'pure' gauge function. Writing 



down the scalar wave equation in the (u, 
solution for (j)(u) as 



v, y, z) system and setting $ = e % 



4>{u), one obtains the general 



ln(0(«)) 



2i 



du 



qg(u) - 27 



ln(qg(u) - 2j) 



duh{u)) 



(3.46) 



For the uniform elec tric f ield, g(u) = Equ + Co with Co being a constant. Both the gauges in Eqn. (3.17) (with 
Eq) and Eqn. ( 3.34 ) (with D\ = D2 = Eq/2) have a similar form for g(u), but differ in the form of h(u). To 



set up a tunnelling scenario the explicit form of h(u) must be known as this will determine whether the system is a 
tunnelling through the barrier or an over the barrier reflection system. For more general forms for g(u), it is clear 
that the gauge invariant quantity E — B ^ which implies non-zero particle production. This is related to the pole 
structure of the first term in Eqn.( 3~46| ) and will be explored further in a future publication. 

In summary, it is seen that by a judicious choice of interpretation of the transmission and reflection coefficients 
in each of the two pure gauges, the standard gauge invariant result can be obtained. This tunnelling interpretation 
works in the case of the mixed gauges too with the gauge parameters deciding whether the system is an 'over the 
barrier reflection' one or a 'tunnelling through the barrier' one. The effective Schrodinger equation that is analysed 
in these gauges is expressed in suitable variables that are combinations of the usual spacetime variables. In cases 
where there is no explicit reduction to an effective Schrodinger form, the exact modes themselves are used to set up a 
tunnelling scenario with the recovery of the tunnelling coefficients. The two mixed gauges that were considered lend 
themselv es to easy solutio ns bu t more complicated ones may be considered in a similar manner. The mode functions 
in Eqn. (3.26) and Eqn. (3.41) are a combination of elementary functions unlike the previously known modes. The 
tunnelling coefficients arise solely due to the presence of a factor s ±lA / 2 in the semi-classical wave functions and in the 
exact mode functions. It is interesting to note that, in the case of black hole spacetimes too, it is precisely a factor of 
the form s ±l£ near the horizon that gives rise to the non-trivial result of Hawking radiation. 

In the last section, we go on to interpret particle production by an electric field in terms of tunnelling between the 
two sectors of the Rindler spacetime. 



IV. COMPLEX PATHS INTERPRETATION OF PARTICLE PRODUCTION IN ELECTRIC FIELD 



In section (III), we noted how particle production can be calculated using the tunnelling interpretation. This 
interpretation gives the same result in both the space as well as the time dependent gauges. The spectrum of particles 
produced by the electric field is not thermal in contrast to the spectrum seen by a Rindler observer. We use the 
formal tunnelling method to show, in a heuristic manner, how this particle production can be obtained by tunnelling 
between the two sectors of the Rindler spacetime. 

The imaginary part of the effective lagrangian ImL e ff for the scalar field in a constant electric field, which is related 
to the probability of the system to remain in the vacuum state for all time, is given by ^ 

j j A 1 (qE ) 2 (-1)»+ 1 / Trm 2 \ 
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where m is the mass, q is the charge and Eq is the magni tude of the electric field. We will derive the above expression 
for ImL e ff using the general arguments given in section (A 1). 

Consider the Hamilton-Jacobi equation for the motion of a particle in an electromagnetic field in (1 + 1) dimensions. 



1 fdF . xt \ 2 1 (OF \ a 1 2 



Qt -2U" # J =° (42) 

where F is the action and — (A*, A x , 0, 0) is the four vector potential. We have neglected the dependence of F 
on the y and z co-ordinates since it does not change the results. In a time-dependent gauge, given by Eqn. ([O]), the 
action F can be easily solved for by using the ansatz F = p x x + f(t) to give, 



F(ti,Xi;to,xo) = Px(xi-xq) ± / dt^J (p x + qE a t) 2 + m 2 (4.3) 



where p x is the momentum of the particle in the x direction. The trajectory of the particle in the (t, x) plane is the 
usual hyperbolic trajectory given by, 

{t -t^-{x- Xi f = -(^j (4.4) 

where ij, Xi are constants. For any fixed position t, there are thus two disjoint trajectories corresponding to motion 
in the two Rindler wedges. 

Let us consider the tunnelling of a particle from one branch of the hyperbola to the other branch and back in the 
imaginary time co-ordinate t (see Fig.l). This means that the particle comes back to the same spacetime point as it 
started from. Choosing the positive sign in Eqn. ( |4.3| ) (this choice give a tunnelling probability that is exponentially 
damped), we have, 



F(t , x ;t , x ) = j> dt^J (p x + qE t) 2 + m 2 = <j> du\/ 1 + u 2 

= i ™L IdrVT^ ~^ = i — r d6 cos 2 (d) 
qE J qE J 

= ^ (4.5) 
qE 

where we have made the following changes of variable (p x + qEot)/m = u = it and r = sin(#). Th e expression for 
exp(z_F) from the above equation is seen to be exactly the same as the exponential term in Eqn. (4.1) for n = 1. The 
same argument can be repeated for the particle tunnelling n times to and fro to give 

to 2 /- 2n7r ,„ iirm 2 , . 

F n {t ,x ;t ,xo) = t—— / d9cos (9) = —^n. (4.6) 
qE J qE 



Again, the quantity exp(iF n ) is seen to match with the exponential part of the nth term in Eqn. (4.1). Therefore, 
the imaginary part of the total effective lagrangian can be written down immediately as 

ImL e ff = (prefactor) exp ( — 7Tm n | (4.7) 
t{ V lEo J 

where the prefactor can only be calculated usin g the exact kernel. However, the dependence of the prefactor on n and 
the phase factor (— 1)™ +1 present in Eqn. (4.1) can be deduced using the following arguments. 

The formal expression of the path integral kernel for the above electric field problem, in the time dependent gauge, 
is given by [g 

K(a,b;s) = (a\e ish \ b) (4.8) 

where K(a,b;s) is the kernel for the particle to propagate between the spacetime points a = {x , x) and b = (y , y) 
in a proper time s and h is the hamiltonian given by 

h = ^{idi - qAi){id l - qA*) - im 2 (4.9) 
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where A 1 is the four vector potential given in Eqn. (3.1) and q and m are the charge and mass of the particle 



respectively. Going over to momentum co-ordinates and considering the coincidence limit x = y, the kernel can be 
written in the form, 



TSt \ * f Q Q 

K(x\y v ;x,x;s) = - / g{x l \y v ;s) 

where Q{x , y°;s) is given by 



(4.10) 



g{x\y°;s) = {x»\e™ H \y ) (4.11) 



and H is the hamiltonian 



H = ~\{w + {Px + qEot)2 + m2 ) =~\{w + q2Elp2 + m ') (4 - 12) 

with p = (t + p x /qEo). In the expression for the kernel in momentum co-ordinates, we have integrated over the 
transverse momentum variables p y and p z . The above hamiltonian is that of an inverted harmonic oscillator. (Since 
all reference to p y and p z have disappeared in H, the dependence of F on the y and z co-ordinates was neglected 
when writing down the expression for the Hamilton- Jacobi equation in Eqn. (4.2).) The expression for the effective 
lagrangian is then given by 

/•OO 7 -1 7>00 7 7>00 

L eS = -i -l K (x ,x ; Xl x; S ) = -——t ™ d Px Q (x° , x° ; s) (4.13) 
Jo s z(zir) j s j_ oc 

We would like to evaluate the propagator Q(x°,x°; s) for a tunnelling situation where the particle tunnels from the 
point x° and back in loops. Since the path integral is not well defined for closed loops, it will have to be evaluated in 

some approximate limiting procedure which is outlined below. 

Since the tunnelling potential is that of a n inverted oscillator, we can use all the results of section ( A 1 b| ) with the 



semi-classical wave functions given in Eqns. ( AlS , A20 ). We would like to first account for the factor (—1)™ that arises 
when a particle tunnels from one side of the barrier to the other and back. Consider an incident wave to the right 
of the barrier and impinging on it. Using the method of complex paths, we rotate this wave in the lower complex 
plane (this is the only route possible for the same reason as that given when rotating a right moving travelling wave 
in the upper complex plane) to obtain a wave again incident on the barrier with a (energy independent) phase factor 
exp(i7r/2) being picked up (other factors dependent on the energy are also picked up but these are not important 
here). Since this wave is moving in the wrong direction, we assume that the particle that has tunnelled through has 
the same amplitude as the rotated wave but is moving away from the barrier. This just involves changing the sign of 
the argument of the exponential in the expression for the rotated wave. Rotating this left moving wave again in the 
upper complex plane now, the final wave obtained is a right moving wave with another extra phase factor of exp(i7r/2) 
being picked up. The total phase change with respect to the incident wave is thus exp(«7r). Since this phase factor is 
independent of the energy, the propagator for the tunnelling process too, after one such rotation, will pick up a phase 
of exp(i7r). Similarly, for n rotations, n taking the values 1, 2, 3, . . ., the phase acquired will be exp(in7r) = (—1)". 
Therefore, the propagator Q for n loops, Q n {x , x°; s), can be written as 

g n (x°, x°; s) = N( Px , m, E)e m7V e lS ^ x0 - x °- s) (4.14) 

where S n (x°,x°; s) is the classical action for the tunnelling process and is the prefactor that arises after evaluating 
the "sum over paths" . This prefactor is not expected to depend on the proper time s since the tunnelling process 
takes place instantaneously or on the number of rotations n. So the only quantities it may depend on are p x , m and 
E. Though the form of N cannot be determined, we can obtain the constraint on N so as to give the correct result 
thereby showing the existence of such a factor. The action for the tunnelling process can be determined by solving 
the Hamilton-Jacobi equation 

1 fdS A 2 1 fas ,„V 1 2 OS n 

2[at + * A ) -2[&- qA ) -2 m+ dl = °- (4 - 15) 

The solution to the above equation is, 

S = -Es +p x x± [ dtyj(p x + qE t) 2 + (m? + 2E) (4.16) 
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Choosing the positive sign and setting p x + qE$t = i^/m 2 + 2_Esin(#), we obtain, 

5 = -J5s + P2; a; ± i / d6>(l + cos(26>)) (4.17) 

2qE J 

For closed paths, with 8 taking the values from to 2nir, the above action can be written as, 

S n = -Es±i m l + J E e (4.18) 

We have thrown away the p x x term while retaining the —Es term since the dependence on the x co-ordinate is really 
irrelevant. Defining a new variable 9 = i9/2qEo and rescaling s = as', one obtains, 

S n = -Eas' + (to 2 + 2E)9. (4.19) 

Choosing a appropriately, S n can be cast into a form that matches the action for a fictitious free particle in (1 + 1) 
dimensions with "energy" aE and "momentum" (to 2 + 2E) satisfying the classical energy-momentum relation, 

aE= ^(m 2 + 2E) 2 (4.20) 

where the particle's "mass" is set to unity for convenience. The above equation determines the quantity a and S n 
can be written in the form 

s „ (M - i;a)= (^ = £fcW (4 . 21) 

where 0\ and O2 are the initial and final states of the free particle with s' being the p roper time taken. (Note that 
(6*2 — &i) = 2imr/2qEo.) Substituting this into the expression for Q n ( x° , x ° ; s) in Eqn. (4.14) and evaluating only the 
integral over s in the expression for the effective lagrangian in Eqn. (4.13) without taking the limits, we obtain 



/ 



A s) = N( Px , to, E)e- [% exp ( to < \~ ^ 

S J & \ £iS. 



= -N( Px ,m,E)e m "—= -— exp {iS n {6 2 A\ *)) • (4.22) 

ia[t)2-9iy 

Notice that the prefactor to th e exponential term has no dependence on the proper time s. Now, we use the form for 
the action given by Eqn. (4. IS) and taking the limits for s from to 00, the effective lagrangian for n loops, L c g(n), 
can be written in the form, 

i {gE f / ttto 2 \ 8 / 2nE \ [°° dp x 

Lcs{n) = 2 (2tt)3 6XP \TW« n ) « 6XP \Wo ) 1-^ ^ ' ] ( } 

where we have set (62 — 61) = 2imr/2qEo. Taking the limit E — > 0, so as to obtain the correct result, and using the 
expression for a in Eqn. (4.20), we find that iV must satisfy the relation 

so that the imaginary part of the effective lagrangian for n loops, L g(n), matches the nth term in Eqn. (|4.l|). 
Therefore, in this manner, the contributions to the imaginary part of the effective lagrangian for the uniform electric 
field can be thought of as arising from the tunnelling of particles between the two Rindler sectors. 



V. CONCLUSIONS 



In conclusion, we see that particle production can be obtained in Schwarchild-like spacetimes in the standard co- 
ordinate systems without requiring the maximally extended manifold. The method of complex paths used in ordinary 
quantum mechanics is modified appropriately to produce a prescription that regularises the singularity in the action 
functional and Hawking radiation is recovered as a consequence. In the case of the electric field, particle production 
in different gauges has been described using the tunnelling description which gives a correspondence between the 
transmission and reflection coefficients and the standard Bogolubov coefficients. The interesting feature of the mixed 
gauges that were considered was that the mode functions could be combinations of elementary functions for certain 
values of the gauge parameters. The method of complex paths also gives a simple interpretation of particle production 
in an electric field as arising due to tunnelling between the two disjoint sectors of the Rindler spacetime. Though we 
have only given a heuristic argument in this paper, we will explore this issue further in a future publication. 
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APPENDIX A: FACETS OF TUNNELLING 



In this section we briefly review the basic concepts of semi-classical quantum mechanics in one dimension and 
formally describe the tunnelling process. We then apply the formalism to two potentials, namely, Vi(x) = —x 2 and 
V%{ x ) = —^/x 2 ! an d calculate the transmission and reflection coefficients for both. 



1. Semi-classical limit of Quantum Mechanics 



Consider a simple one dimensional quantum mechanical system with an arbitrary potential V(x) where x denotes 
the space variable (see Ref. [pi for details). To describe the transition of the system from one state to another, we first 
solve the corresponding classical equations of motion and determine the path of transition. This path is, in general, 
complex since many processes like tunneling through a potential barrier cannot occur classically. Therefore, the 
transition point qo where the system formally makes the transition is a complex number determined by the classical 
conservation laws. Then, the action S for the transition from some initial state Xi to a final state Xf given by 

S(xf,x l ) = S{x f ,q ) + S(q ,x l ) (Al) 

is calculated. Here, S(qo,Xi) is the action for the system to move from the initial state Xi to the transition point qo 
while S(xf, qo) is that to move from go to Xf . The probability P for the transition to occur is given by the formula 

P~exp (-^Im[S{x f ,q )+S(q ,Xi)]^ . (A2) 

The above formula is valid only when the argument of the exponential is large. Further, if the potential energy has 
singular points, these must also be considered as possible valu es f or qo. If the position of the transition point is not 



unique, then it must be chosen so that the exponent in Eqn. (A2) has the smallest absolute value but still must be 
large enough so that the above formula be valid. 

If the transition point qo is real, but lies in the classically inaccessible region, then the above formula gives the 
transmission coefficient for penetration through a potential barrier, while if the transition point is complex, it solves 
the problem finding the over the barrier reflection coefficient. The ~ sign in the above formula is used since the 
coefficient in front of the exponential is not determined. This can be determined by calculating the exact semi- 
classical wave functions. Generally, it is desirable to find the ratios of two different transitions so that this coefficient 
does not matter. 

The physics of the tunnelling and the "over the barrier" reflection processes are very different. In the tunnelling 
process, the semi-classical analysis gives a transmission coefficient that is an exponentially small quantity with the 
corresponding reflection coefficient being unity. In contrast, in the "over the barrier" reflection process, just the reverse 
is obtained. The transmission coefficient is unity while the reflection coefficient is an exponentially small quantity. 
Both these processes will be encountered when the electric field is studied in different gauges. 

We will now review the method of calculating the transmission and reflection coefficients for a typical quantum 
mechanical problem using the method of complex paths for a general potential V(x). 



a. Description of the method of complex paths 

Consider the motion of a particle of mass to in a region characterised by the presence of a potential V(x) in one 
space dimension. The problem is to calculate the transmission and reflection coefficients between two asymptotic 
regions labelled L and R where the semi-classical approximation to the exact wave function is valid. After identifying 
these regions and writing down the semi-classical wave functions, definite boundary conditions are imposed. The usual 
boundary conditions considered are that in one region, say L, the wavefunction is a superposition of an incident wave 
and a reflected wave while in the second region R, the wave function is just a transmitted wave. Then, a complex path 
(in the plane of the now complex variable x) is identified from R to L such that (a) all along the path the semiclassical 
ansatz is valid and (b) the reflected wave is exponentially greater than the incident wave at least in the latter part of 
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the path near the region L. The transmitted wave is then moved along the path to obtain the reflected wave and thus 
the amplitude of reflection is identified in terms of the transmission amplitude. Having done this, the normalisation 
condition is used i. e. the sum of the modulus square of the transmission and reflection amplitudes should equal unity, 
to determine the exact values of the transmission and reflection coefficients. 

For a given potential, the turning points go (or transition points) are given by solving the equation 



p(q ) = y/2m(E-V(q Q )) = 



(A3) 



where p(x) is the classical momentum of the particle and E is the energy of the particle. In general, go is complex. 
At these points, the semi-classical ansatz is not valid since the momentum is zero. Further, the potential can possess 
singularities. At these points too, the semi-classical approximation is invalid. Therefore the contour between the two 
regions should be chosen to be far away from such points. In general the contour will enclose them. Therefore, the 
relation between the transmission and reflection amplitudes is determined by taking into account the turning points 
and the singularities of the potential. 

The Schrodinger equation to determine the wave function ip of the particle is 



2m dx 2 



(E-V(x))ip. 



(A4) 



Referring to ||, the semi-classical wave function, in the classically allowed region where E > V(x), is given by the 
formula 



ip = Cip 1/2 exp ( - / p(x)dx j + C 2 p 1/2 exp . 



p(x)dx 



(A5) 



while in the classically inaccessible regions of space where E < V(x), the function p{x) is purely imaginary and the 
wave function is now given by the relation 



<l> = C x \p\ 



-1/2 



exp 



- / \p(r)\d.r 



C 2 \p\ 



-1/2 



exp Q J \p(x)\dx 



(A6) 



The condition on the potential for semi-classicality of the wave function to be valid is 



dx \p{x) 



or, in another form, 



mh\F\ 



<1, 



< 1, 



F = — 



dV 
dx 



(A7) 



(A8) 



It ought to be noted that the accuracy of the semi-classical approximation is not such as to allow the superposition 
of exponentiall y sm all terms over exponentially large ones. Therefore, it is inapplicable in general to retain both 
terms in Eqn. (|A6| ). We will consider a few cases of interest in this paper and refer the reader to Q for a exhaustive 
discussion along with suitable illustrative examples. 



Consider the case in which the semi-classical condition (A8) holds in the regions x — ► ±oo. As x — > — oo, the wave 
function is assumed to be a superposition of incident and reflected waves and is written in the form, 



ip=p 



-1/2 



exp 



p{x)dx) +C 2 p- 1/2 



exp 



p(x)dx 



where the incident wave has unit amplitude while the reflected wave has amplitude given by C 2 . As x 
wave function is assumed to be a right moving travelling wave, 

ip — -Jz exp ^ — J pdx 



(A9) 
-oo, the 

(A10) 



The method of complex paths can now be applied on the function (AlO). The contour is chosen either in the upper 
or lower complex plane such that the reflected wave is always exponentially greater than the incident wave along that 
part of the path near the region x — ► — oo. If this is satisfied along one of the contours then C 2 is determined in terms 
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of C3. To carry out the above procedure however, the exact semi-classical wave functions as x — > ±00 have to be 
determined. This will be done explicitly later for the relevant cases. 

A different case arises when the semi-classical ansatz holds in the vicinity of the origin x — rather than at x = ±00. 
The boundary conditions assumed in this case are the same as above with the condition x — > 00 replaced by x — > + 
and x — ► —00 by x — ► _ . Here, the required contour is about the origin and is chosen to be small. But it must still 
be large enough so that the reflected wave is much larger than the incident wave along the latter part of the contour 
near the region x < 0. 

In the above cases the method of complex paths gives the exact transmission and reflection amplitudes. But, 
in certain cases it is enough to assume that the transmission and incident amplitudes are equal to unity while the 
reflection amplitude is exponentialy damped and consequently very small. Here, the "over the barrier" reflection 
coefficient for energies large enough so that the reflection coefficient is exponentially small, has to be determined. 
In this case, the condition E > V{x) is always satisfied. Therefore, the transition point go at which the particle 
reverses its direction is the complex root of the equation V(qo) = E. Let go lie in the upper complex plane for 
definitencss. Now, the amplitudes of the incident and transmitted waves are equal (both are set to unity within 
exponential accuracy). To calculate the reflection coefficient, the relation between the wave functions far to the right 
of the barrier and far to the left of the barrier must be determined. The transmitted wave can be written in the form 



I I 1 

exp 



Vp 



pdx 



(All) 



where x\ is any point on the real axis. We follow the variation of ipx along a path C in the upper complex plane which 
encloses the turning point go. The latter part of this path must lie far enough to the left of qo so that the error in the 
semi-classical incident wave is less than the required small reflected wave. Passage around go only causes a change in 
the sign of the root \J~E — V(x) and after returning to the real axis, the function ipx becomes the reflected wave ipR. 
Going around a complex path in the lower complex plane converts into the incident wave. Since the amplitudes 
of the incident and transmitted waves may be regarded as equal, the required reflection coefficient is given by 



R 



exp 



Im / pdx 



(AI2) 



Now we can deform the contour in any way provided it still encloses the point go- In particular, the contour can be 
deformed to go from x% to go and back. This gives 



R = exp 



4 

-— Im 

n 



pdx 



(AI3) 



Since p(x) is real everywhere, the choice of X\ on the real axis is immaterial. The above formula determines the above 
the barrier reflection coefficient. It must be emphasised that to apply the above formula the exponent must be large 
so that 1 — R is very nearly equal to unity. 

Finally consider a situation where the amplitudes of the reflection and incident wave are equal. The transmission 
coefficient is now an exponentially small quantity. This case corresponds to the standard tunnelling process. The 
problem is characterised by the presence of real turning points between which lies the classically forbidden region 
where the energy E < V(x). For definiteness, let there be two real turning points labelled g_ and q + . The potential 
V(x), in the immediate vicinity of the turning points of the barrier, is assumed to be of the form 



E-V(x)*F (x-q±), 



F 



dV_ 

dx 



(A14) 



This assumption is equivalent to saying that the particle, near the turning points, moves in a homogeneous field. 
With this assumption, the amplitude of transmission C3 is given by (refer to Ref. || page 181) 



C 3 = exp I - 



9+ 



\p(x)\dx 



(AI5) 



The transmission coefficient is then given by | C3 1 2 . The above formula holds only when the exponent is large. In 
the derivation above, we have assumed that the semi-classical condition holds across the entire barrier except in the 
immediate vicinity of the turning points. In general, however, the semi-classical condition need not hold over the 
entire extent of the barrier. The potential, for example, could drop steeply enough so that Eqn. ( A14) is not valid. In 
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these cases, the exact semi-classical equations have t o be determined before applying the method of complex paths. 
The cases encountered in this paper all satisfy Eqn. ( A14 ). 

We now apply the above results to two potentials. The first is the well known inverted harmonic oscillator potential 
V\(x) — —g\x 2 with g\ > while the other is V2(x) = —g^jx 2 with g 2 > 0. The first potential arises when the 
propagation of a scalar field in a constant electric field background is studied. The second potential arises when the 
propagation of a scalar field in Schwarzchild-likc spacetimes is considered in the vicinity of the horizon. 



6. Application to the potential V\{x) 



Consider the potential given by 



Vx{x) 



-gi% 



(Ai6) 



where g\ > is a constant. This potential is the inverted harmonic oscillator potential and is discussed extensively in 
many places (see for example Rj , || ) . We will follow the semi-classical treatment given in Ref. || and briefly review 
the calculation of the reflection and transmission coefficients for both the tunneling and over the barrier reflection 
cases. 

The semi-classicality condition ( |A8| ) for the above potential is, 



hgix 



/ 2m[Ei + gxx 2 } 



213/2 



< 1 



(A17) 



where m is the mass of the particle and E\ is its energy. The above condition definitely holds for large enough \x\ 
and for any value of E±, either positive or negative. Therefore the motion of a particle moving under such a potential 
is semi-classical for large enough and hence holds as x — > ±00. 

Since the motion is semi-classical for large \x\, we can expand the momentum p(x) as, 



p(x) = \j2m {Ei + gix 2 ) « y / 2mg 1 (x + — — 

\ 2gn 



(A18) 



Using Eqn. (AlS), the semi-classical wave functions can be written as follows with the following boundary conditions. 
As x — > 00, we assume that the wave function is a right moving travelling wave ipn while as x — * — 00, it is a 
superposition of an incident wave of unit amplitude and a reflected wave given by tp^. Therefore, we have, 



ipL - (-£)~ l£l ~^e- l « 2/2 + C 2 (-0' lei '^e^ 2/2 
where we have made the definitions 



-00 ) 
-00) 



2m<7i 



1/4 



£1 




Ei 



(A19) 
(A20) 



(A21) 



Following the variation of Eqn. ( A19| ) around a semi-circle of large radius p in the upper half plane of the now complex 
variable £, we obtain 



C 2 = —iCa cxp (— 7rei) 
The conservation of particles is expressed by the condition that, 



|C 3 | 2 + |C 2 | 2 = 1 



From Eqn. ( |A22| ) and Eqn. ( A23 ), the transmission coefficient is, 

1 1 



T = C 3 
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(A22) 
(A23) 

(A24) 



while the reflection coefficient is 
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R = \C 2 ? 



(A25) 



Note that the passage through the lower half complex plane to determine C 2 is unsuitable since on the pa rt of the 
path —ir < <j> < — ir/2 , where cf> is the phase of the complex variable £, the incident wave (first term in Eqn. ( A20 )) is 
exponentially large compared with the reflected wave. The above formula h olds for a ll energies E\. This is because, 
even for negative energies, the semi-classical wave functions given in Eqns. (A19 A20) are exactly the same with the 
boundary conditions being fully satisfied. 

If Ei is large a nd ne gative, Eqn. flA24|) give s T w e -^V 2m /9i\ E i\/ n anc j thus R ~ 1. This is in accordance with the 
formula in Eqn. (A15). To apply Eqn. (A15) it is necessary to calculate the turning points first. The real turning 
points are q = — y/\Ei\/gi and qi 



V\ E i\/9i- Therefore, 
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This gives the same answer. 

We can calculate the over the barrier reflection coefficient using Eqn. (A13) for E\ large and 'positive. The turning 
points now are given by solving the equation p(qo) — with the condition that E\ > V\{x) always. Since E\ > 0, the 
turning points are go = 'E\/ g\. Choosing the positive sign for qo and setting x\ = 0, the integral in Eqn. ( |A13| ) 
is evaluated as follows: 



p(x)dx = \j2mgi j y/E\jg\ + 



'la 



/■! l 

i\/2mjg\E\ I y/l - y 2 = -iir V '2m j ' g\E\ 
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(A27) 



Therefore, 



R — exp ( — — iry / 2m/giEi 



(A28) 



The above formula can also be obtained directly from Eqn. (A25) by neglecting the exponential term compared to 
unity which means that the energy has to be large. 



c. Application to the potential Viix) 
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Consider the potential given by 
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(A29) 



where g 2 is a positive constant. The potential has a singularity at the origin. This potential arises when the effective 
Schrodinger equation is calculat ed for Schwarzchild-like spacetimes. This aspect will be dealt with in later sections. 
The semi-classical condition (A8) for this potential takes the form 



%2 



< 1 



(A30) 



where E 2 is the energy. It is clear that the above relation holds for large |x|. It also holds for small \x\ if \J2mg 2 3> h. 
Notice that the quasi-classicality condition for small |x| is independent of the sign and magnitude of the energy £'2- 
For this potential, we will be concerned only with the small |x| behaviour in contrast with the potential V\. Since the 
motion is semi-classical for small |x|, we expand the momentum p(x) as 
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p(x) = x j2m[E 2 + ^|) « v + xli^Ezx (A31) 



Notice the similarity between Eqn. (|Alq ) and Eqn. ( A31 ). 

We will calculate the over the b arrier reflection coefficient with E 2 > and small, which will be of interest later. 
Using the expansion in Eqn. ( A31 ), the semi-classical wave functions with the following boundary conditions, namely, 
that for x > the wave function is a right moving travelling wave while it is a superposition of an incident wave of 
unit amplitude and reflected wave for x < 0, are 

i, R = C 3 r 2+i e 4 « 2/2 K > 0) (A32) 

V^l = (-£)-'^ + i e -^ 2 /2 + C 2 (-0 l£2+i e l Z 2/2 (£ < 0) (A33) 

where we have made the definitions 



Following the variation of Eqn. ( A32) around an small semi-circle of radius p < \jg2J\E2\ ( m contrast to the potential 
V\ where the radius p was large) in the upper half complex plane, we obtain, 

C 2 = C 3 exp('-7r£ 2 + y) (A35) 

Setting T — |C 3 | 2 = 1 and R = |C 2 | 2 , we finally obtain, 

R = Te- 27rE2 = r e -s27rVS^ (A36) 
Using the normalisation condition R + T = 1, we obtain, 

1 e -i-27rV2mg2 

T = J == and R = 7 == (A37) 

Notice that the above result is independent of the energy E% and hence holds for E2 < too. For small the lack 
of dependence on Ei is not too surprising since the contour is such that it is not too close to the real turning points 
go = ± y/ g2 1 1 Ei I . Anyway, when E2 ~ + , p is "large" and therefore the contour is chosen to lie in the upper complex 
plane for the same reason as given in the a nalysis of the pote ntial V\ in the previous section. 

We will derive the result in Eqn. ( |A36| ) using Eqn. ( A12 ). The complex turning points qo are the roots of the 
equation E2 = —92/^0 where E 2 > and therefore, the turning points are go = H\/g2/ E 2 — ±ipo- Hence, we have 
to evaluate the integral, 



J pdx = v / 2mE 2 J y 1 + ^dx (A38) 

where the contour C encircles the point x — ipo in the upper half complex plane. H owev er, since there is a singularity 
at the origin, we cannot deform the contour as was done when deriving Eqn. (A13). Therefore, as a means of 
regularisation, we modify the potential to 

V mod (x) = (A39) 

X z + € £ 

where the limit e — ► must be taken at the end of the calculation. The turning points for the modified potential are 
a^mod = ±z-\/e 2 + g 2 /E 2 while the poles of the modified potential are at x — He < x mo d- Even in this case, there is 
a singularity on the path of integration which contributes to the integral rather than the turning point. Therefore, 
integrating upto +ie using the modified potential and back, we obtain, 



J pdx = lim 2^2mE 2 J yl 

= lim 2i \J2mE2t I dy 
Jo 
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o V 1 - y 
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= iit \j 2mE 2 po = iir\j2mg 2 (A40) 
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We therefore recover the result given in Eqn. (A36). From the above calculation it is clear that, due to the singularity 
at the origin, the reflection coefficient has no contribution from the turning point at all. 
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FIG. 1. AB and CD are the hyperbolic branches representing motion in the two disjoint Rindler sectors for a fixed t. The 
particle tunnells from B to C and back in a circular path along the imaginary time-space plane. If the particle makes n-loops 
around the circle, the imaginary part of the action will be n times the value for single loop. 
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